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Abstract 

We apply the recently developed method of differential renormalization to the Wess- 
Zumino model. From the explicit calculation of a finite, renormalized effective action, 
the /^-function is computed to three loops and is found to agree with previous existing 
results. As a further, nontrivial check of the method, the Callan-Symanzik equations 
are also verified to that loop order. Finally, we argue that differential renormalization 
presents advantages over other superspace renormalization methods, in that it avoids 
both the ambiguities inherent to supersymmetric regularization by dimensional reduc- 
tion (SRDR), and the complications of virtually all other supersymmetric regulators. 
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1. Introduction - Since it has been developed, differential renormalization has had 
some important successes. At considerable calculational (and conceptual) simplifica- 
tion over other regularization methods, it has been used to fully renormalize massless 
\ip 4 theory to three loopsfl], to calculate the triangle anomaly [1,2], and the 1-loop 
/3-function in Yang- Mills theory [1]. Further, mass effects in \ip A theory have also been 
included straightforwardly [1,2], and the translation of the method to a setting in which 
counterterms appear explicitly has also been achieved in the course of proving the 
unitarity of the method for massless \<p theory [3]. Like for the \(p 4 model, the per- 
turbation theory of the Wess-Zumino model is fairly simple and well-known, but at the 
same time sufficiently nontrivial for it to represent a good testing ground for differential 
renormalization, and we find that the method works as well in this case as it does in 
\ip A theory. 

The basic idea of the method is as follows: rather than regularizing by altering 
the field content or the propagators of a theory (e.g., Pauli-Villars, higher derivatives, 
point splitting) or the space on which it is defined (e.g., dimensional regularization), 
it regularizes and renormalizes only the x-space amplitudes that are too divergent to 
have a Fourier transform into momentum space. In practice, this is done by writing 
amplitudes in x-space, expressing the divergent pieces as derivatives of less singular 
terms, and then formally performing partial integrations freely (which in fact will be 
precisely the implicit regularization and renormalization that the method performs). 
Counterterms do not appear at any point (although it is possible to keep track of 
them explicitly [3]), and one goes directly from bare amplitudes to renormalized ones, 
without a separate intermediate step of subtraction in the regularized theory. The sim- 
plest example occurs at one loop both for the 4-point function in Xip A theory and the 
2-point function in the Wess-Zumino model: these amplitudes will contain the term 
((x — x') 2 ) -2 , and its Fourier transform into momentum space (of course, in the end 
one must go to momentum space to calculate scattering amplitudes) is: 

r e'i'P-x 4tt 2 r°° 1 
J x 4 p J H 

where r = (x 2 )2 ; and J\ is the first-order Bessel function. For r — > oo the integrand 
converges fast enough, but for r — > 0, J\(pr) ~ \pr, and we have a logarithmic UV 
divergence. To regularize it, differential renormalization then prescribes the use of the 
following identity: 
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valid everywhere except at the origin. Now, the term on which the D'Alembertian is 
acting is, as opposed to a well-defined distribution with a proper Fourier trans- 
form. By formally integrating by parts we are then able to define a regularized Fourier 
transform for \: 



J x 4 A J x 2 A J 
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where M = and 7^ = 1.78107... is the Euler-Mascheroni constant. A general fea- 
ture of differential renormalization which can already be seen from the example above 
is that in the process of expressing singular functions as derivatives of less singular func- 
tions and then integrating by parts, thereby discarding (infinite) surface terms, we are, 
on the one hand, removing the divergence of an amplitude and, on the other, acquiring 
a new mass parameter M which appears in the guise of an integration constant. Since 
the method does not keep track explicitly of the infinities being subtracted, this mass 
parameter will be crucial for the study of the RG invariance of renormalized amplitudes 
and the consistency of the method. 

Because Eq.(l) is an identity except at the origin, bare and renormalized ampli- 
tudes will also be identical except at isolated points. The method does not, so to speak, 
uniformly and blindly regularize an entire theory ab initio (like the regularization pro- 
cedures mentioned above), but rather removes divergences singly where they appear in 
amplitudes. In this sense, differential renormalization is a "minimal" regularization: it 
alters the theory in the least amount possible, and in particular it does not disturb its 
original symmetries. 

To apply it to the Wess-Zumino model, we will use the manifestly supersymmetric 
Feynman rules of Grisaru, Rocek, and Siegel[4,5] (except for the fact that we do not go 
into momentum space), to calculate loop corrections to the chiral propagator. For su- 
perspace theories in general, using these Feynman rules, one ends up with the following 
superspace effective action: 



r= f ateY. I d 4 x 1 ...d 4 x N F 1 (x 1 ,0)F 2 (x 2 ,9)...F N (x N ,e)G^ N \x 1 ,...,x N ), (2) 
J N J 

where the Fi(xi,9) are functions of superfields and their (supersymmetry) covariant 
derivatives, and are translationally invariant functions of i spacetime coordinates. 
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Because of the non-renormalization theorem [4, 5], the Wess-Zumino model, in particu- 
lar, does not have a genuine coupling constant renormalization (but only the one inher- 
ited from wavefunction renormalization) [6,7], and is entirely renormalized by removing 
the divergences from GW(x - x'). That we do, up to three loops, with differential 
renormalization; with more ease than previous methods, we find a three-loop /3-function 
which agrees with the existing results[7,8,9], and as further consistency checks of the 

(2) 

method, we verify the Callan-Symanzik equation for Gf-en, at each loop order, also to 
three loops. 

2. The Model and Loop Calculations - The Wess-Zumino model is described by the 
following superspace action: 



S = Jd 4 x (| " | (/ + J d W 3 )) , (3) 



where (x a , 9 a ,0 a ), a = 1, 2, 3, 4, a = +, — , a = +, — are coordinates of d = 4, N = 1 
superspace, and <fi (<fi) is a chiral (antichiral) superfield: 

In component form, 

(j)(x, 9, 0) = <p(x) + Oil>{x) + 66F(x) + i6o- a 6d a ip(x) - % -eed a ip{x)o a 6 + ^OOOOUipix) 



and 

s = J A 



-d a <p*dP<p - % -{$o a d a ^ - d a ^) + FF* - + X -^A + c.c.) 



For our calculations in superspace, we use the Feynman rules of Grisaru, Rocek and 
Siegel[4,5]; in coordinate space, they are: %) the propagator is 

<^(x,9)cP(x',9')> = 1 5^(6-6') 

4:ir z (x — x'Y 



Att 2 (x - X 1 ) 2 



it) for each vertex, include a factor of 4; Hi) for a chiral (0 3 ) vertex, include two factors 



3! 

of —\D 2 acting on any two of the propagators for the lines arriving at the vertex; for 
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vertices containing an external line, include only one — \D 2 factor acting on any of the 
two other (internal) lines; iv) same for antichiral vertices (<^> 3 ), with — \D 2 factors; v) 
integrate over internal superspace points (xi n t,0i n t) and external 0's, 9 ex t] vi) include 
appropriate symmetry factors for each diagram. Finally, the D-algebra and transfer 
rules are: 



D 2 D 2 D 2 = 160D 2 , D 2 D 2 D 2 = 16QZ) 2 , 

= D 3 = , [D Q , D 2 } = 4ia a ^d a DP , 



D a {e',d x ,)[8^{e-9')f{x-x')} 



-D a (9,d x )[5^(9-9>)f(x-x')}. 



With these rules for calculating Feynman diagrams, the effective action takes the 
form (after all ^-integrations but the last have been carried out): 



d?6 



d^xdy^x, 9, 9)4>(x', 9, 9)G^ (x — x 



+ 

+ 
where 



J d 2 9 J d 4 x 1 d%d%<P(x 1 ,9,0)0(x 2 ,9,9)0(x 3 ,9,9)G < $(x 1 ,X2,x 3 ) + c.c. 



d^9 



(4) 



G^ 2 \ 



x — x 



S^\x — x') + K\{x — x') + k 2 {x — x') + k 3 (x — x') + etc. (5) 



G^(xi,x 2 ,x 3 ) = S^Hxi - x 2 )5 ( - A \xi - x 3 ) + bperms. 
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(6) 



and the dots in Eq.(4) refer to higher-point functions which are either finite or are 
renormalized by the renormalization of G^ 2 \ The diagrams contributing to G^ are 
depicted in Figs.(l) and (2). n\ is the one-loop contribution from Fig.(l.a), k 2 the two- 
loop contribution, Fig.(l.b), and K3 contains contributions from the four three- loop 
diagrams in Fig. (2). 

The non-renormalization theorem determines that the renormalization of g comes 

(3) 

solely from wavefunction renormalization, and dictates the above form for G cf . L , correct 

(3) 

to all loops. By studying the Callan-Symanzik equation for G ch , we find a simple 
relaton between the /3-function and the anomalous dimension 7(g) of <p: 



= (p{g)^-^rig)jG^(xi,X2,x 3 ;g) = 
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=> (3{g) = 3g 7 (g) . (7) 
The Callan-Symanzik equation for the renormalized 2-point function then becomes: 



M Wl + m Tg ~ 2l{3) ) G "" {X ~ X ' ! S ' M) = 

( M w + "w(|-^)) G ~ (x - x ' ;9 ' M) - °- <8) 

In general, M is a renormalization mass scale that governs the RG flow of the the- 
ory; in differential renormalization this scale appears through the process described in 
the introduction. 

We now proceed to renormalize the theory. For the diagrams of Figs.(l) and (2) 
the .D-algebra and ^-integrations are straightforward, if tedious, and we present only 
the final x-space expressions to be renormalized: 

Ki{x - x) = jxhr (9) 

"2(i^) (x - X 1 ) 2 l ^ (X - Xl ) 2 ( Xl ~ X'Y (10) 



K 2 (X - X 1 ) 



1 ( 9 \ 6 1 



/4 f 4 1 1 1 1 

d xi J d X2 . _ , A , _ s 2 , _ , 2 , _ ^,, 2 



K 3c (x - x') 



K 3d (x - X') 



x\) 4 (xi - X2) 2 (x - x 2 ) 2 (x 2 - x') 2 

1 x 
4 \Att 2 ) {x-x') 2 

I ^ I ^ (X - Xl ) 2 ( Xl - X 2 )H X 2- X'Y (13) 



1 ( 9 V 1 

2 Uvr 2 / (x-x') 2 

r d 4 Xi f d 4 X2 1 1 1 1 1 

J J (x — xi) 2 [x' — X2) 2 [x\ — X2) 2 (x' — xi) 2 (x — X2) 2 



X 



This last (finite) integral was done with the aid of the Gegenbauer polynomial 
technique [1,7]. Now, in the course of performing the above integrals, we use the 
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following differential renormalization identities [1] whenever singular terms appear (i.e., 
terms like the l.h.s. below): 



1 lnx 2 M 2 
■-□ ^ 



lnx 2 M 2 
In 2 x 2 M 2 



1 In 2 x 2 M 2 + 2\nx 2 M 2 
8 x z 

1 ln 3 x 2 M 2 + 31n 2 x 2 M 2 + 61nx 2 M 2 
— □ ! 

12 x 2 



x^O. 



(15) 
(16) 
(17) 



With this, we are able to calculate the renormalized values of Kf. 



K\ en {x 
K T 2 en {x 

kT(x 
K r 3 ?(x 



2 / 1 \ 4 lnx 2 M 2 

4 /l\ 6 ln 2 x 2 M 2 + 21nx 2 M 2 
9 UJ ° x 2 



6 \A-k 



^ ' 47T ' =i 



□ 


X 2 




1 
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In 3 x 2 M 2 + 2 In 2 x 2 M 2 


+ 4lnx 2 M 2 




x 2 




□ 


ln 3 x 2 M 2 +3 In 2 x 2 M 2 


+ 61nx 2 M 2 



3 9 V4vr. 

> . , fi / 1 \ 8 lnx 2 M 2 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 



where we have set x' = for simplicity. To calculate the /3-function and verify the 
consistency of the method we apply the Callan-Symanzik equation, Eq.(4), to G^en = 

5 (4) + K ren + ft ren + K ren_ Wg write 

/%) = Pig 3 + fog 5 + fog 7 



and separate the pieces in the Callan-Symanzik equation according to order in g 2 and 
the coefficients of the 5-function and the different powers of log. All of these have to 
vanish separately. At 0(g 2 ), we get the one-loop /3-function: 
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At 0(g 4 ), the coefficient of S(x) vanishes for the following value of the two-loop 
/3-function: 



(25) 



and the coefficient of □ ln X J^ vanishes for the above value of (3\; this is a consistency 
check of the method. At 0(g 6 ), the vanishing of the coefficient of 5{x) gives the following 
three-loop /3-function: 



/tf= (| +9((3) ) (^)V. 



(26) 



Finally, the coefficients of □ ln and D ln x ^ m2 will vanish for the above values 
of /?i and /?2, and these are then two further consistency checks. Our final result for 
the /3-function is then: 





( 9' 








V47T, 




V4vr) 



4 21 
+ (^+9C(3)) 



4 
4vry 



(27) 



We note that although the three-loop part is scheme-dependent, as opposed to 
the one- and two-loop results, the coefficient of £(3) should be universal (barring, of 
course, a coupling constant redefinition involving £(3)) because the only diagram at 
three loops leading to the transcendental function £(3) is primitively divergent. While 
our coefficient of £(3) coincides with that obtained in [9], it does not agree with [7]. The 
correct result is that of [9], however [10]; in fact, because that is a four loop calculation, 
it is possible to verify the three-loop result by RG pole equations (which the authors 
of [9] in fact do). 

Conclusions - We have seen that the method of differential renormalization is read- 
ily applicable to the superspace Wess-Zumino model and, in particular, the /3-function 
calculations were considerably simpler than in SRDR [7,9]. Firstly, the integrations 
performed here were either trivial (integrations with <5-functions) or very simple con- 
volutions; secondly, we did not have to keep track of subtraction of sub divergences: at 
each loop order the renormalization is done in a single step, by the direct use of "dif- 
ferential renormalization identities" like Eqs.(12)-(14) whenever singular expressions 
are encountered (in SRDR, on the other hand (cf. Eq.(2.11) of [7]), a careful account 
needs to be kept of different renormalization constants corresponding to contributions 
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of sub divergences at higher loops). 

We have attempted to describe the differential renormalization of the Wess-Zumino 
model in such a way that the extension to other superspace theories becomes obvi- 
ous: one should again start with the general form for the effective action, Eq.(2), with 
the GW written in terms of x-space integrals, and apply differential renormalization 
identities to these integrals, thus regularizing and renormalizing the theory. Differ- 
ential renormalization will take translation-invariant quantities G^ into translation- 
invariant quantities G^2, and will manifestly maintain super symmetry. Naturally, we 
expect this procedure to be far simpler than several other existing regularization meth- 
ods (e.g., the supersymmetric versions of Pauli-Villars, higher derivative and point- 
splitting methods). Furthermore, if the method proves to be as simple to implement 
in general as SRDR, it would have the advantage that the ambiguities associated with 
SRDR would be avoided: like for dimensional regularization, SRDR presents ambigu- 
ities due to the presence of intrinsically four-dimensional quantities, like e^ vpc , in a 
dimensionally continued setting [5,11]. Differential renormalization, of course, never 
leaves four dimensions and would thus avoid these problems altogether. One initial 
indication of the further applicability of the method is the fact that gauge (vector) 
superfield propagators are identical to the massless chiral propagators we have consid- 
ered here. Thus, calculations with vector superfields should not in principle present 
any new difficulty; for instance, the one-loop correction to the vector propagator in 
supersymmetric Yang-Mills theory coupled to a massless chiral superfield is again as 
simple to compute as in SRDR. Further work along these lines is in progress. 

Acknowledgments - It is a pleasure to thank J.I. Latorre for always illuminating 
discussions. Comments on the manuscript from D.Z. Freedman and J.I. Latorre are 
also very appreciated. 

References 

[1] Freedman, D.Z., K. Johnson and J.I. Latorre, to be published in Nucl.Phys. B. 
[2] Freedman, D.Z., Proceedings of the Stony Brook Conference on Strings and 

Symmetries, Spring 1991. 
[3] Freedman, D.Z., R. Muhoz-Tapia, X. Vilasfs-Cardona, manuscript in preparation. 
[4] Grisaru, M.T., M. Rocek and W. Siegel, Nucl.Phys. B159(1979)429. 
[5] Gates, J., M.T. Grisaru, M. Rocek and W. Siegel, Superspace, Benjamin-Cummings, 
1983. 

[6] Iliopoulos, J. and B. Zumino, Nucl.Phys. B76(1974)310. 

Ferrara, S., J. Iliopoulos and B. Zumino, Nucl.Phys. B77(1974)413. 



9 



[7] Abbott, L.F. and M.T. Grisaru, Nucl.Phys. B169(1980)415. 

[8] Townsend, P.K. and P. van Nieuwenhuizen, Phys.Rev. D20(1979)1832. 

[9] Sen, A. and M.K. Sundaresan, Phys.Lett. B101(1981)61. 

[10] M.T. Grisaru, private communication. 

[11] Siegel, W., Phys.Lett. B94(1980)37. 

Figures 



Fig.(l): One- and two-loop contributions to G( 2 \x — x') 



Fig.(2): Three-loop contributions to G ( - 2 \x — x'). 
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